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We present a numerical investigation of the dynamics of symmetry breaking in both Abelian and 
non-Abelian [577(2)] Higgs models in three spatial dimensions. We find a class of time-dependent, 
long-lived nonperturbative field configurations within the range of parameters corresponding to 
type-1 superconductors, that is, with vector masses (m„) larger than scalar masses (m s ). We argue 
that these emergent nontopological configurations are related to oscillons found previously in other 
contexts. For the Abelian-Higgs model, our lattice implementation allows us to map the range of 
parameter space - the values of = (m a /m v ) 2 - where such configurations exist and to follow them 
for times t ~ O(10 5 )m _1 . An investigation of their properties for z-symmetric models reveals an 
enormously rich structure of resonances and mode-mode oscillations reminiscent of excited atomic 
states. For the SU (2) case, we present preliminary results indicating the presence of similar oscillonic 
configurations. 



I. INTRODUCTION 

The mechanism of spontaneous breaking of local gauge 
symmetries plays a fundamental role in our current un- 
derstanding of high-energy particle physics [1] and of con- 
densed matter systems [2]. In very general terms, it can 
be stated that a spontaneously broken symmetry is al- 
ways associated with the existence of degenerate vacuum 
states: the theory predicts the existence of discrete or 
continuously degenerate vacua while Nature chooses one 
of them. Strictly speaking, the mixed matrix elements 
describing possible transitions between vacuum states 
only vanish in the infinite- volume limit. However, since 
mixed matrix elements scale with volume V as exp[— cV], 
where c is a positive constant, for all practical purposes 
symmetries do get broken for large enough volumes. 

In the realm of relativistic quantum field theories, the 
particle spectrum of a given model is computed as small 
perturbations about a broken-symmetric vacuum state. 
In addition to these, models with nontrivial nonlinear 
couplings may also have nonperturbative, solitonic so- 
lutions to their equations of motion [3]. These usually 
come in two kinds: for models with nontrivial vacuum 
structure, there can exist static solutions that owe their 
stability to the topology of the vacuum manifold, the 
so-called topological defects [4, 5]. For example, models 
with a Id (one spatial dimension) real scalar field with a 
double- well potential have kink solutions, while Abelian- 
Higgs models in 2d have Nielsen-Olesen vortices. In both 
cases, symmetry is restored at the core of the topological 
defect. In contrast, nontopological solitons owe their sta- 
bility to the conservation of a global charge [6, 7]. The 
distinctive signature of both topological and nontopolog- 
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ical solitons is their time-independence: they are effec- 
tively static solutions to the equations of motion, even 
in the case of Q-balls, where the complex scalar field 
is written as $ ~ exp[— iu)t] so as to transform the time- 
dependent term in the equation of motion to a mass term 

Given the vast richness of spatiotemporal phenomena 
in Nature [8] , one should suspect that other possible non- 
perturbative configurations exist in relativistic field the- 
ories once we allow for their time-dependence. Usually, 
these are not taken into account, as they are expected 
to be short-lived and hence dynamically uninteresting. 
However, if long-lived configurations do exist, they are 
bound to play a crucial role when fields are far from equi- 
librium and, in particular, during symmetry breaking. If 
present, they will change our understanding of the vac- 
uum, as they comprise, together with possible topological 
and nontopological extended field configurations (EFCs) , 
nonperturbative fluctuations about it: for example, if we 
are to sum over possible contributions to the path inte- 
gral, these must be included. In a cosmological setting, 
if we are to study the approach to thermalization during 
post-inflationary reheating [9, 10], these configurations 
may possibly be very important [11, 12]. 

As more recent research has shown, there are abun- 
dant examples of such long-lived solutions. [13, 14, 15, 
16, 17, 18, 19, 20]. The first hint was the discovery 
of breathers in Id kink-antikink scattering (see, e.g., 
ref. [21]): for certain relative velocities a new time- 
dependent configuration, a breather, would form. Re- 
markably, breathers were never seen to decay. Their 
demise may only occur through highly-suppressed non- 
perturbative decay modes. Even before that, spherically- 
symmetric configurations called "pulsons" were found by 
Bogolyubsky and Makhankov [22, 23]. These were the 
first examples of the configurations later called oscillons 
in ref. [13, 14], where it became clear that such oscilating 
large-amplitude real scalar field solutions are present in 
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any 3c? model with amplitude-dependent nonlinearities as 
long as the potential has a region of negative concavity, 
V" < 0. For more details on the properties of scalar field 
oscillons see ref. [24]. 

The more recent extension of oscillons to the Stan- 
dard Model presented in refs. [15] and [18] prompted 
us to search for oscillons in the Abelian-Higgs model. 
In our first work [19], we found that vortex-antivortex 
annihilation in 2d can, for a range of parameters, gen- 
erate remarkably long-lived oscillon-like configurations. 
These are characterized by a see-saw oscillation in the 
^-component of the magnetic field and seem to owe their 
stability to a gauge-field induced mass gap for the scalar 
field, although an in-depth study is still lacking. 

In the present work we explore the existence of oscillon- 
like states in 3c? Abelian-Higgs models. We find that they 
not only exist but are quite easily found in the context 
of type-1 superconductors, that is, for gauge-field masses 
substantially larger than Higgs-field masses. As in 2d, we 
construct an effective phase diagram mapping the range 
of parameters where oscillons are found. Using extrap- 
olation, we obtain an approximate critical value of the 
control parameter (3 = (m s /m v ) 2 beyond which we con- 
jecture that no oscillons are produced during symmetry 
breaking. 

The crucial point, though, is that these configurations 
emerge spontaneously during symmetry breaking. We 
don't start with an approximate spherically-symmetric 
solution and see it relax into oscillons, as has been the 
rule in such studies. The U(l) oscillons literally condense 
dynamically, as the system transitions from its symmet- 
ric to its broken-symmetric state. Although it is true 
that we take advantage of the formation of flux-tubes to 
facilitate the formation of oscillons - a 3d analog of the 
vortex-antivortex annihilation in 2d - those flux tubes 
occur naturally during the symmetry-breaking process. 
It is quite remarkable that, in a model where there is 
no topologically-stablc defect, flux-antiflux tubes annihi- 
lation will form long-lived EFCs. Even though they do 
not show the same symmetry restoration at the core as 
do topological defects, we show that in 3d U(l) oscillons 
clearly probe into the V" < part of the potential. We 
also briefly discuss the remarkably rich structure of these 
objects, which display resonant mode-mode fluctuations 
reminiscent of excited atomic states. Finally, we present 
preliminary results in the context of SU(2) Higgs models 
where, although we were not able to find long-lived con- 
figurations, there is a strong indication that they exist in 
the type-1 regime. 

This paper is organized as follows: in section II we 
present the model and our conventions. In section III 
we describe our search for U(l) oscillons and how we 
managed to isolate them to investigate their longevity. 
In section IV we present the phase diagram for 3d U(l) 
oscillons and obtain the critical value of the control pa- 
rameter (3. In section V we describe the rich resonant 
structure of these configurations, stressing the similari- 
ties with excited atomic states. In section VI we present 



our preliminary results for type-1 SU(2) models. In sec- 
tion VII we present our conclusions and a brief summary 
of our results. Finally, the two appendices describe the 
technical details and issues of the lattice implementation 
of Abelian and non-Abelian models. Particular attention 
is paid to the proper implementation of gauge constraints 
in the presence of stochastic forcing terms. 

II. 17(1) & SU(2) EQUATIONS AND 
CONVENTIONS 

We use the continuum Lagrangian, 

C = X> M *t • X>"* - -F" v ■ + -A($ f • $ - T] 2 ) 2 , (1) 

where 2? M = d^ + igA^, and is either the Abelian U(l) 
or non-Abelian SU(2) gauge field. The context within 
which we will be using one or the other will be clear. 
The operation • does nothing for U(l) but is equal to 
a-b= iTr a)b for SU (2) matrices. Performing the scaling 
A^ — ► r]~ 1 A fl , cf> — * rf 1 ^ and x —> "qgx, there is only 
one independent parameter in these models, the ratio 
of scalar to vector masses, (3 = (m s /m v ) 2 = X/(2g 2 ). 
For convenience, we will keep A = 2 and vary the gauge 
coupling g so that all times are quoted in units of the 
scalar mass, [t] = to" 1 . 

Note that this choice of conventions can be related to 
condensed matter models that use the hopping parameter 
k, (for example, see reference [25]) by K = l/y/2g. The 
critical parameters separating type-1 (small (3) and type- 
2 (large (3) for {7(1) superconductivity is /? = 1 and K — 
1/V2- 

The continuum equations for this system are, 

= d&V; (2) 

V* = J v , (3) 

where J M = i (^D^ - $D>*&) is the conserved cur- 
rent. We will solve these equations numerically in the 
temporal gauge Aq = 0, which makes the time compo- 
nent of the last equation a nondynamical constraint. This 
gauge also allows the definition of simple functions for the 
conjugate momenta and hence for a Hamiltonian lattice 
implementation. Details of the lattice implementation 
can be found in the Appendices. 

III. FINDING OSCILLONS IN THE U(l) 
THEORY 

In a previous work, we showed that in 2c? it is possible 
to find long-lived, time-dependent oscillon configurations 
from the annihilation of vortex-antivortex pairs (hence- 
forth vav) [19]. These remarkable EFCs are characterized 
by a persistent see-saw oscillation of the magnetic field B z 
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FIG. 1: (Color online) A few observables characterizing the 
2d (7(1) oscillon. The blue line is the maximum energy density 
of the EFC while the black line is the total energy within a 
radius of R = 4. The lower plot shows both the expected 
radius of the EFC R e ff (green) and the minimum value of 
$*$ within the EFC (red). All of these observables except 
the energy are oscillatory in time and have a nontrivial fourier 
space distribution of frequencies. 



and very little emitted radiation. [The interested reader 
can see an animation here.] We were also able to show 
that the formation of U(l) oscillons can be described as 
a phase transition in field configuration space with (3 as 
the control parameter: oscillons were shown to form from 
vav annihilation only for (3 < (3 C ~ 0.13(6) ± 2. We con- 
structed an order parameter E osc /E v , where E osc is the 
oscillon energy and E v ~ 2tt(3 1 ^ is the Nielsen-Olesen 
vortex energy, showing that E osc /E v ~ (|/3 — {3 C \)°' 2 ^ ±2 . 

Figure 1 shows a few relevant observables during the 
oscillon phase. Note the constancy of the total energy 
(black line), integrated over a finite volume surround- 
ing the oscillon. The blue line denotes the local max- 
imum of the total Hamiltonian density, max.[H(t,x,y)]. 
The program scans the lattice at each time step to find 
the local maximum of the energy density. The dashed 
green curve displays the value of the effective radius of 

the configuration, computed as R 2 f f = A. r„ J^}\ > f r , with 

Jo H(r)dr 

rn ~ few mj 1 . The red curve denotes the minimum am- 
plitude of the scalar field, mm[^^(t,x,y)]. Here, it is 
important to notice that the vacuum is at $t$| vac = 1, 
while the inflection point is at $^$|i n f = 0.5. Thus, the 
2d oscillon is marginally nonperturbative as the field am- 
plitude hovers just around the inflection point. We will 
see that this will not be the case in 3d. 

It is natural to search for similar oscillon-likc EFCs in 
3d. However, as soon as this is attempted, one meets 
a few obvious challenges. First, contrary to 2d, there 
are no stable topological vortices in 3d. We can find 



(7(1) flux tubes, but they do not share the stability of 
2c? Nielsen-Olesen vortices. Second, it is computationally 
much harder to search for EFCs in 3d. As we will see, if 
we start from a symmetric initial state and quench to the 
broken symmetric state, there will be an excess of energy 
that makes potential candidates harder to isolate. 

Generally, finding an oscillon involves some method of 
confining enough energy within a small spatial region and 
then seeing if the energy remains localized without any 
or with very little dissipation to infinity. This method 
implies that oscillons can be thought of as attractors 
in field configuration space: through their natural dy- 
namical evolution, the interactions among the fields will 
conspire to create an oscillon solution. This property of 
oscillons has recently been demonstrated analytically in 
the context of pure real scalar field models [24]. A simi- 
lar proof for gauge models is lacking, but our numerical 
results indicate that this property will also hold for these 
more realistic models. 

From our experience with the 2c? theory we know that 
for small j3 = m 2 / m 2 , oscillons can form from vav decay. 
We can try to extend this result to 3ci. Note that the 
fact that oscillons form easily from vav annihilation for 
small p doesn't mean that they only form for small (3: 
we have observed that they can also form from a zero- 
phase Gaussian initial condition in <f> for larger [3, $ = 
1 — $o exp[— r 2 /R 2 ), where $o is the departure from the 
vacuum at |<I>| 2 = 1 and R its spatial extension. 

We stress that there is a fundamental difference in 
these two approaches. Using a Gaussian as an initial con- 
dition relies on our knowledge of the approximate oscillon 
behavior and is equivalent to coaxing the solution into ex- 
istence by letting the fields relax into it from a nearby 
point in field configuration space. This is commonly done 
in numerical relaxation techniques. A quench, on the 
other hand, doesn't assume any initial profile for the so- 
lutions: the oscillons emerge spontaneously as the sys- 
tem works to minimize its energy and maximize its en- 
tropy dynamically. Thus, finding oscillons through this 
dynamical approach offers strong evidence of their exis- 
tence as attractors in field configuration space. It means 
we should expect them to be present during symmetry 
breaking. 

We now describe the procedure to form oscillons from 
vav annihilation. We first thermalize the fields in a sym- 
metric single well potential at low temperature T and 
with viscosity 7 = 1. We do this using the stochas- 
tic Langevin approach described in the Appendices. We 
then switch to the double well potential, while at the 
same time turning off the stochastic noise but maintain- 
ing the viscosity. The viscosity will dampen excessive 
energy and will allow a more transparent identification 
of the many vortices and anti- vortices formed on the lat- 
tice. In steps: 

• Following the lattice implementation in appendices 
1 and 2, set all fields (A^,®) and their derivatives 
to zero. 



4 



• Thcrmalize using Langevin dynamics in a quadratic 
potential with minimum at <f> = 0. The tempera- 
ture T should be chosen to generate a distribution 
of fluctuations across fc-space. A typical number is 
T ~ 0.1m s . Using 7=1 (in units of m s ), the sys- 
tem should thcrmalize well within t ~ 10m J 1 . We 
take the system to be thermalized when d t (II 2 ) ~ 0, 
where II = dt® and the brackets denote a volume 
average. 

• Switch to T — > while keeping 7=1 and simulta- 
neously switch to the double well (Higgs) potential. 

• Then evolve for t ~ 15m -1 until the vortices and 
antivortices form. 

• Set 7 = 0, and then evolve conservatively. 

Nearby vortices (anti-vortices) will then interact to 
form higher-TV vortices (anti-vortices) while nearby vav 
pairs annihilate. In 2d, for strongly type-1 parameters, 
(j3 < 0.136) the process of vav annihilation has a good 
probability of forming gauged oscillons [19]. Whether 
they form this nontopological, time-dependent, radially- 
symmetric bound state depends on both the size of the 
perturbative modes around the structure and on the rela- 
tive velocity of the vortices as they interact. Our strategy 
in 2d was to minimize as much as possible their relative 
velocity using the viscosity. We also noticed the robust- 
ness of the formed oscillons against perturbative radiative 
modes. 



A. Attempt 1 in 3d: Gaussian ansatz 

In 3d the dynamics is significantly different, as the flux- 
tubes are not topologically stable as are vortices in 2d. 
Also, there are more decay paths, and more surface en- 
ergy to compensate for in higher dimensions. We will 
attempt both the Gaussian ansatz initial condition and 
the stochastic field quench as we did in 2d to search for 
oscillon-like structures. As an aside, we note that our 
approach to find new time-dependent nonperturbative 
solutions is applicable to a variety of nonlinear partial 
differential equations. The same way that oscillons could 
never have been predicted with the perturbative analyt- 
ical techniques commonly used to study nonlinear equa- 
tions, we conjecture that many nonperturbative time- 
dependent solutions remain unknown. Small-amplitude, 
spatially-extended oscillons were recently found in the 
context of real scalar field theories a posteriori, after they 
were numerically discovered. These objects are amenable 
to a small amplitude treatment, as has been shown in 
refs. [11] and [26]. 

In 3d, the first result we report is a negative one. In 
figure 2 we plot the total energy (blue, upper curve) and 
its contribution from gauge fields (red, lower curve) for 
a Gaussian perturbation from the vacuum. The ansatz 

is $ = 1 — (S> en* , with <f> = 1 and R = 4 and for 
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FIG. 2: (Color online) Attempt to find a 3d U(l) oscillon from 
a Gaussian initial condition for f3 = 0.04 (or g = 5). Plotted 
is total energy (blue, upper curve) and magnetic field energy 
(red, lower curve) within a 10 3 box. Note that, although the 
configuration is unstable, the leaking of energy is much slower 
than expected in a linear theory. 



(3 = 0.04. No information was put in the gauge fields, 
that is, we took Aj = 0. We tried many different com- 
binations of radii, amplitudes and couplings g, but none 
of the configurations were stable. Even though this same 
type of ansatz was successful in the equivalent 2d the- 
ory for a large range of (3, it is clearly too simplistic for 
3d. As we know that 3c? U(l) oscillons exist (see be- 
low), a longer-lived oscillon should be found with the 
ansatz method, although one would need to incorporate 
the gauge fields in a more sophisticated way. As we are 
more interested in the dynamical emergence of oscillons, 
we will not pursue this further. It is worth mentioning 
that although the configuration decays in t < TOOmJ 1 , 
this is still enough time for the object to dramatically 
affect macroscopic physics. It is important to keep this 
in mind when judging the possible implications of short 
lived nonperturbative resonances. 



B. Attempt 2 in 3d: Breaking the Symmetry 

Our next attempt involved making a string network 
and letting those strings interact to see if long-lived, lo- 
calized nontopological spatiotemporally-complex objects 
were formed. To do this, we followed the same steps 
outlined above: first, we thermalized the field to Ti att = 
T8x~ d in a quadratic potential; second, we quenched it 
by switching to the double-well potential to seed the for- 
mation of strings; third, we maintained the viscosity long 
enough so that the strings stabilized and the scalar field 
approached its vacuum expectation value. The viscosity 
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dynamically decreased the temperature T — > 0. As in 
2ci, the only consideration when choosing an initial T is 
to generate enough excited modes to seed the formation 
of a sufficient density of strings upon symmetry break- 
ing while not so high that they are overdense. Specific 
parameters are given below. Finally, after the string net- 
work formed, the friction was turned off and the strings 
interacted. We then had to search for nontrivial struc- 
tures resulting from these interactions. We started to ex- 
plore at very small 0, since in 2d that is when structures 
are more likely to form. Also, we expect f3 c to be smaller 
than in 2d, as there will be more surface tension to com- 
pensate for. In what follows, we describe the details of 
two searches, characterized by the initial parameters as 
configurations 1 and 2, respectively. This should allow 
our results to be reproducible by other groups. 

Configuration Ci 

The details of configuration C\ were: 

{d, N L , L,Sx, 5t, g, X,T],T latt } Ci = (4) 
{3,96,19.2,0.2,0.05,4,2,1,0.25}, 

where the various symbols stand for: d - spatial di- 
mensionality; Nl - number of lattice points; L- lattice 
length; Sx - lattice spacing; St - time discretization; g - 
gauge coupling; A - scalar coupling; r\ - viscosity; Ti at t 
- lattice temperature. As remarked, the basic procedure 
for constructing this configuration is similar to how we 
formed vortices in 2d, although we have more sensitive 
dependence on the time-scales for string formation in 3d. 
For clarity, the steps are: 

• Set all fields to their vacuum configurations, $ ~ 
in the quadratic potential V = <1>^<I>. 

• Thcrmalizc using Langevin dynamics with 7=1 
for t ~ lOra^ 1 at r latt = 0.25. 

• Set T — (no Langevin kicks), 7 = 0.25, and si- 
multaneously switch to the double well (Higgs) po- 
tential. Evolve for t = SmJ 1 . 

• Evolve with 7 = 1.0 for t ~ 12to s _1 until string 
loops form. 

• Switch 7 — > 0, and satisfy condition in eq. A13 
(the Gauss constraint) by setting all the momenta 
to zero. 

• The remainder of the run is now energy-conserving 
and we evolve this portion with an 0(St s ) symplec- 
tic integration scheme. 

Configuration C\ nicely generates an initial string net- 
work. As shown in figure 3, due to the small lattice size, 
the strings annihilate quickly into a loop which then vi- 
olently decays into a few oscillon candidates. Because 
some of these candidates are relativistically boosted away 



from the string annihilation region, they are also Lorentz- 
contracted in the direction of motion. Some of these 
decay and we are left with a fairly hot background 
with three spatially-localized configurations show in fig- 
ure 4. One of these objects still moves quite fast and 
is noticeably Lorentz-contracted (top left of figure 4) ; it 
decays after colliding with the larger slow-moving one 
which has quickly taken a spherical shape (center right 
of figure 4). The second fast structure then decays at 
t ~ 200mJ 1 (not visible in figure), after interactions with 
the hot background, while the spheroidal object remains 
until t ~ 1000m" 1 . This is our oscillon candidate. 

There are a few points to address. First, we call the 
fast-moving, localized object (on the top left of the fig- 
ure 4) an oscillon candidate and not a radiation wake 
because we observed that it does not decay according 
to an expected dispersion relation in the direction of 
travel or perpendicular to it. One would expect a wake 
to disperse at least perpendicularly to the direction of 
travel. The full simulation shows that one can eas- 
ily distinguish and identify the radiation wakes. View 
simulation here. 

Second, if these relativistically-boosted structures de- 
cay in about ~ 100 — 200to~ 1 time units, does it mean 
they are not very stable? It is hard to say in the context 
of this simulation. We have seen in 2d that even config- 
urations that will evolve to become oscillons, forming at 
rest and without a hot background (in a near vacuum), 
can take a long time O(10 2 mJ 1 ) to settle down into a 
coherent object. Thus, if they were boosted with relative 
velocities to a large amplitude radiation bath at birth or 
soon after, they may never have settled into an oscillon 
state. Of course, since this is a Lorentz invariant theory, 
we can boost any configuration across the vacuum with- 
out affecting its stability. But if the boost exposes the 
object to a relatively fast moving and large amplitude set 
of radiation modes (so, not near- vacuum) , its decay may 
be catalyzed. Here is an intuitive justification for this 
fact. Consider a configuration-space attractor character- 
ized by a coherent object. Consider further that it is pos- 
sible to attach a measure in configuration space, so that 
there is a well-defined "distance" from any configuration 
to this attractor point. In the presence of large perturba- 
tions, the probability of an initial configuration settling 
into the attractor will decrease with its "distance" to it. 
Back to our simulation, we see that candidate configu- 
rations far from reaching the oscillon state may not be 
stable enough to settle into it in the presence of large per- 
turbations. However, our surviving configuration clearly 
does. 

The third point is that any stability arguments have 
to take into account stochastic thermal effects. This par- 
ticular configuration (radiative fields over whole volume 
plus EFC contribution) has an average energy-density 
(7i) ~ 0.06267 and an approximate temperature of the 
same order of magnitude. While the approximate energy 
density of the oscillon is more than an order of magnitude 
higher in the core, much of the lower-energy modes which 
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FIG. 3: (Color online) 4 snapshots (t = {12, 18, 24, 43} m a 1 increasing left to right) for configuration Ci. Plotted is energy 
density with an isosurface of Hi = 0.275 (blue), $ f $ = 0.45 (red), \{E 2 + B 2 ) = 0.275 ( green). (Energy isosurface is 25% 
transparent and when outside of red, appears purple). At the first time-slice the friction has just been turned off. As the string 
network twists and accelerates, it generates violent changes in the magnetic field at its sharp corners. The shrinking of the 
toroid forces the flux to get radiated out of these sharp edges, generating the magnetic wakes seen in the last time slice. View 
simulation here. 



make up the structure will be disrupted by such a tem- 
perature. As has been shown in the context of real scalar 
field oscillons, it is not surprising that thermal noise will 
compromise the lifetime of coherent states [27]. 

Even though no similar stability analysis has been per- 
formed in the context of gauged oscillons, temperature 
effects may be a serious obstacle to the formation of os- 
cillons and to their stability. This is not really that sur- 
prising, since we expect symmetries to be restored at high 
temperatures. The crucial question, then, is if the gauged 
oscillons can sustain stochastic thermal effects. In other 
words, if they are present, at what fraction of T c do they 



get destroyed? 

Back to our simulation, after ~ 250771J 1 only one os- 
cillon is left on the lattice. It can be treated as a sin- 
gle object with some residual velocity in the midst of a 
thermal bath. To characterize this object we can look 
at a few observables. The maximum energy-density is 
max[7i(i, x, y, z)] ~ 4. The fraction carried by the gauge 
fields is quite large at some times, H^ifgs ~ 3. Simi- 
larly to the 2d gauged oscillons, it is clear that the gauge 
fields play a fundamental role in the oscillon dynamics. 
The maximum amplitude of scalar field oscillations is 
$t$ - 0.27. This is very important: although sym- 
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FIG. 4: (Color online) Residual oscillon candidates from con- 
figuration Ci after the radiation from the flux tube decay has 
dissipated throughout the lattice. Time slice at t = 74mJ 1 . 
Energy density isosurface is at 7i± = 0.275 (blue); scalar field 
amplitude is at = 0.45 (red not visible since it is in- 

side the energy isosurface); gauge fields energy density is at 
\{E 2 + B 2 ) — 0.275 (green). 

metry is not quite restored at the oscillon core as it is 
in usual topological defects, it still qualifies as a large 
nonperturbative coherent fluctuation away from the vac- 
uum state, since the inflection point of the potential is at 
^^linf = 0.5. This particular configuration then dies at 
t ~ lOOQm-V 

Configuration Ci 

The details of configuration Ci are: 

{d, N L ,L,6x,St,g, A, r?, Ti att } Cl = (5) 
{3,96,19.2,0.2,0.05,5,2,1,0.25}. 

We prepare this simulation very similarly to configura- 
tion C\. the gauge coupling is increased from g = 4 
to g — 5 and we also extend the time in which the 
dissipation is present for formation of string loop from 
t ~ \2rn~ 1 to t ~ 26m7 1 . 

• Set all fields to their vacuum configurations in the 
quadratic potential <&^<i>. 

• Thcrmalize using Langevin dynamics with 7=1 
for t ~ IOto" 1 at ILtt = 0.25. 

• Switch T — > and 7 — > 0.25 and simultaneously 
switch to the double well (Higgs) potential. Evolve 
for t = 3TTIJ 1 . 



• Evolve with 7 = 1.0 for t ~ 26m s 1 until string 
loops form. 

• Set 7 = 0, and satisfy condition A13 by setting all 
the momenta to zero 

• The remainder of the run is now energy-conserving 
and we evolve this portion with an 0(5t 8 ) symplec- 
tic integration scheme. 

An oscillon is formed from this configuration and we 
observed it for a time of t ~ ISOOOmJ 1 . At this point we 
stopped the simulation as the oscillon had not changed 
for thousands of time units, continuing to drift across 
the lattice. In figure 5 we plot for a short time the evolu- 
tion of two important observables related to this oscillon: 
maximum of energy-density max[7i] and minimum am- 
plitude of scalar field min[<I>^<I>]. Compare with similar 
plot for equivalent oscillon in 2c? in figure 1 . Cleaner data 
will be obtained below. 

Now that we know that these objects exist in 3d, we 
can move on to a more in-depth study of their properties. 
The problem with the two previous configurations is that 
they leave too much energy on the lattice, making it hard 
to study the oscillon accurately. This extra energy also 
catalyzes their decay. Next, we propose an approach to 
remove all the spurious energy so that we can investigate 
(7(1) oscillons in more detail. 

C. Isolating the 3d U(l) Solution 

In order to isolate the oscillon, we take an initial con- 
figuration that we know allows for at least one from flux- 
tube annihilation, that is, for the annihilation of two long 
strings with opposite magnetic fluxes inside, the equiva- 
lent of 2d vortex-antivortex annihilation. After the flux 
tubes interact and annihilate, we search the lattice for the 
maximum value of the energy density, max[7i(t, x, y, z)] 
which hopefully correlates with the presence of an oscil- 
lon. After finding the maximum and tracking it for a 
while to make sure it is sufficiently long-lived, we place 
a spherical friction wall at a radius r > Rf from it (the 
choice of Rf to be made explicit soon) to dissipate all 
modes outside the EFC, effectively "mopping" the lat- 
tice from possible destructive perturbations that may af- 
fect the formation and longevity of the oscillon. The 
fields within the cavity at r < Rf settle into an oscil- 
lon, which moves about the lattice carrying the friction 
wall along. When the energy for r > Rf is only very 
slowly decreasing we know that we have dissipated most 
of the excess energy in the lattice volume; any residual 
energy within the cavity is contributed by the oscillon it- 
self. Once the average energy density outside the cavity 
has decreased to ~ O(10 -4 ) times the oscillon's energy 
(that is, the energy for r < Rf), we can turn the friction 
off and watch the oscillon drift about the lattice. Using 
this method, we have seen oscillons live for longer than 
t > 69000mJ 1 for — 0.04 without showing any signs 
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FIG. 5: (Color online) On the left is plotted the minimum value 
the right is the maximum of the energy density TL(t, x, y, z). 

of instability. In figure 6 we show the oscillon energy 
as a function of time for a simulation with (3 = 0.01, 
where its long time behavior can be explicitly seen. Note 
that even though our phase diagram stops at (3 = 0.04 
this doesn't mean that there are no oscillons formed for 
smaller (3. This why we include the example in the fig- 
ure with (3 = 0.01. The small (3 analysis is incomplete 
due to CPU limitations and not physics, as is the case 
for (3 > f3 c . This method also allows the approximate 
determination of the oscillon's energy, obtained by a vol- 
ume integral over the lattice [E osc (f3 = 0.04) ~ 9], and 
the value of the energy at the oscillon's core, which is 
the tracked value max[7i(i, x, y, z)]. This value and the 
value of the minimum of <i>^<i> are plotted in figure 5 for a 
short time after the oscillon has formed and the friction 
wall had been turned off (marking the t = in the plots.) 
These plots correspond to the isosurface plots in figure 7. 

More specifically, the friction is implemented by con- 
structing a shell around the location of maxpi. (t, x, y, z)] 
which vanishes for a radius r < Rf — 4.0 and equals 
7 = tanh(r — -R/)/2 for r > Rf. The oscillon is safely lo- 
calized within this shell. Although the implementation of 
the friction wall does a great job isolating the oscillon so- 
lution, it obviously violates the Gaussian constraint A13, 
as one can see by rederiving current conservation from 
the equations of motion with this nontrivial form of 7. 
However, the fact that what violates Gauss's constraint 
is a viscosity term only acting in a region with negligi- 
ble energy-density contribution from the configuration of 
interest is helpful. If the violations are far enough from 
the region of interest - that is, for r < Rf - then the 
simulation is trustworthy since the violations continue to 
decrease in magnitude as r increases. We also note that 
usually, when there are numerical instabilities or viola- 
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FIG. 6: (Color online) Time evolution of the total energy 
(dashed) and of the energy within a sphere of R = 4 for a 
f3 — 0.01 oscillon. After a quick shedding of a large amount 
of energy, the onset of the oscillon stage at at time t ~ 10 3 is 
characterized by a very slow energy decay accurately fitted by 
E(t) — —c±t~ C2 , where C2 — 1.2. The inset shows the numeri- 
cal value of the energy and the fitted curve (indistinguishable) 
for t > 10 4 . 



tions in constraint equations, these will tend to create 
exponential instabilities. This is not the case here, as 
our numerical method is perfectly stable in the presence 
of these nonpropagating charge densities. 

To investigate the violation in Gauss's constraint we 
take the absolute value of equation A13. This we then 
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FIG. 7: (Color online) 4 snapshots (t = {0, 1, 16, 90} m" 1 ) of the initial stages of a long-lived 3d U(l) oscillon. Plotted is total 
energy density (blue) at an isosurface of 0.05. Energy density isosurface in the gauge fields (green) is also at 0.05, while for the 
scalar field magnitude <& = 0.9 (red). Size of plotted area is 1906mJ 3 and f3 = 0.04. The energy of the resulting oscillon is 
about E osc ~ 9. The spherical shell (magenta) in lower plots denotes the friction wall which follows the oscillon throughout the 
lattice. Because of the toroidal boundary conditions, the nearly spherical profile of the low energy outer surface of the oscillon 
is visible as a cross section at t = 90 (lower right). Note also that a second oscillon candidate is present at t — 1 (top right), 
but doesn't survive long due to the friction. View simulation here. 



plot along with the friction wall and the localized struc- 
ture which has formed inside of it, the oscillon. Hopefully, 
from this experiment we can gain some understanding of 
what effect, if any, the Gaussian constraint violation has 
on the dynamics. In figure 8 we show the oscillon (in 
blue an green) up to t ~ dOOmJ 1 as it moves through re- 
gions of Gaussian violation (in orange) at about the 10 -3 
level without any discernible variation of its observable 
max[7i] and min|<i>^<i>| oscillations as compared to oscil- 
lons formed in systems with Gaussian violations at the 
10~ 12 level of absolute magnitude. As the friction wall is 



turned off at t ~ 2QQm~ 1 (bottom left), the extra charge 
density from the Gaussian violation stays fixes in space, 
while the oscillon keeps moving about the lattice with 
no discernible effect (bottom right). View simulation 
here . 

In figure 9 we plot the expectation value 
of the violation of Gauss's constraint A(t) = 
(\h^M-Ei_ l )-2gIm(ni<Z>)\) as as function 

of simulation time. At 200m J 1 we turn off the spherical 
wall of friction and then any violation left on the lattice 
is frozen due to conservation of charge to our numerical 
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FIG. 8: (Color online) 4 snapshots it — {1, 45, 192, 520} mj 1 ) of 3d U(l) oscillon displaying violation of Gaussian constraint. 
Plotted is energy density (blue) with an isosurface of 1.00 marked for scale. Energy density in the gauge fields (green) is also 
marked at 0.10 while <&' <& = 0.4 (red). The absolute value of the violation of the Gaussian constraint is shown in orange at an 
isosurface of 5 x 10 -3 . Size of plotted area is 1906m J 3 . This is also a /3 = 0.04 configuration. The energy of the resulting oscillon 
is about E oac ~ 9. The spherical shell (purple, bottom plots) denotes the friction wall that follows the oscillon throughout 
the lattice. After t > 2007T17 1 we turn the friction wall off, but plot its last position to highlight that the region of Gaussian 
violation stays anchored in space. View simulation here. 



precision. 



Based on these observations, we did not see any un- 
desirable effect due to the presence of small amounts of 
violations in A(x, t). What we gain from this relaxation 
of numerical rigor is a path to investigate the EFC for a 
long time. 



IV. CONSTRUCTING A 3d 17(1) PHASE 
DIAGRAM IN CONFIGURATION SPACE 



Once we have numerically shown the existence of a 
time-dependent EFC in a given theory, the next step is 
to categorize the region of parameter space in which this 
solution is an attractor configuration. From our expe- 
rience, it is almost impossible in some regions to actu- 
ally find an object, despite hints that it might exist. So, 
the method we chose to predict where the object lives 
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FIG. 9: (Color online) Global expectation of the violation 
= (\^J2 i (Ei-E i ._ l )-2qIm(Ui^)\) for the Simula- 
tion in figure 8. For simulations without a spatiotemporally 
complex support function for friction, A(t) ~ W~ 12 (at) b , 
where a and b are small numbers. 



in configuration space is to pick an observable which can 
be used to characterize certain properties of EFCs in a 
region where they are known to exist, and then study 
that observable as a function of parameter space. With 
this we can predict by extrapolation the regions with no 
EFCs. As in 2d, the results can be usefully organized in 
a phase diagram. 

The observables we choose must vanish when no 
oscillons are present. Since, when oscillons are 
present, the maximum local energy-density on the lat- 
tice max[7i(i, x, y, z)\ is much larger than the background 
noise, and because we are searching for large amplitude 
structures, we will use the maximum energy-density and 
the minimum min[<I>^<I>] amplitude as order parameters. 

Our basic method will be to form an oscillon follow- 
ing similar steps to configuration Ci in section IIIB, ex- 
cept that we will repeat the procedure for various fi. On 
the lattice, we track max[7i] and min[<I>T<I>] as a function 
of time, and take a running average of these oscillating 
observables. As they converge to a near-constant value 
(usually at times t > lC^mJ 1 ), we plot them in a phase 
diagram as shown in figure 10. We can then use the di- 
agram to extrapolate and find the critical parameter (3 C 
beyond which oscillons do not form. The basic method 
can be summarized as follows: 

• Prepare a configuration similarly to Ci in section 
IIIB. If an EFC is formed (confirmed by identify- 
ing a persistent local maximum in energy density) , 
then measure max[7i] and min[<i>^<f>] to confirm it is 
an oscillon. If not, repeat with a different Langevin 
realization using the same initial temperature. 

• If an oscillon forms then introduce the moving fric- 



FIG. 10: (Color online) Phase diagram for the 3d U(l) 
flux-antiflux tube annihilation — > oscillon transitions. Plot- 
ted as a function of /3 are the order parameters Oe(/3) = 
i?max(/3)/-Emax(0) with numerical values in boxes, and 
O s (/3) = (1 - $ t $ min (/3))/(l - $ + $ mi n(0)) in circles. The 
continuous lines are fits which we use to extrapolate to the 
critical value where no oscillons are expected to form by this 
mechanism. We find p c ~ 0.0893 and p c ~ 0.0908 for the 
Oe{0) and 0$(/3) order parameters, respectively. 



tion wall until the energy-density outside the co- 
herent object is negligible compared to max[7i]. 
(t ~ 250mJ 1 has proven to be sufficient.) 

• Turn off the friction and allow the system to evolve 
using an 0(St s ) symplectic integration routine. 

• Track max[H] and min^*] for 2500m;r 1 . 

• Construct a running average of these quantities 
calling them -E max and <&T(j> min respectively. If this 
running average converges to a near-constant value, 
then take this value as a point in the phase diagram. 

It is essential that within a reasonable amount of time 
the time-averaged observables show only negligible time 
variation. If we cannot (after some reasonable amount 
of attempts) find an oscillon, we move to a lower /3 since 
our intuition from the 2d work indicates this favors their 
production. Fortunately, there is a regime in which the 
observables converge well in a short time, allowing us to 
fit a function and extrapolate to the critical point in (3 
above which no oscillon should form. 

In figure 10, we present the averaged order parameters 



o B (P) 



__ ffmax(/3) 
^max(O) 



(6) 



where E max ((3) (blue curve) is the value take from the 
running average of max[7i], and in the green curve we 
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FIG. 11: (Color online) On the left is plotted the minimum value of for an j3 — 0.04 oscillon configuration. The right is 
the maximum energy density TL. These correspond to t = 2 x W^mJ 1 in data and highlight the very slow leaking of energy 
from oscillons. When min^O] > l/s/2 we expect the oscillon to decay completely. The red lines are 3-parameter fits to the 
data which are identified in equations 8 and 9. 
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The numerical values are fitted from the region where we 
trust the data and then extrapolated to an effective P c . 
We find /3f ~ 0.0893 and pf ~ 0.0908 for the energy and 
$ order parameters, respectively. These are higher than 
the region in which we actually see the oscillons form. 

Take, for example, f3f obtained from the fit to $t$. 
We extrapolated to when the min[<I>^<I>] — + 0, thus match- 
ing the vacuum. This is not restrictive enough. Gener- 
ally, these oscillons decay if their core oscillations do not 
probe near the inflection point of the potential. If we 
had extrapolated to the inflection point $ ~ l/\/2, Pc 
would change to (3 C ~ 0.078. This more restrictive cri- 
terion moves f} c to almost exactly where we stop seeing 
oscillons from flux tube decay. 

In figure 1 1 , the time evolution of these (unsealed) ob- 
servables is presented for j3 — 0.04. This configuration 
is at about the limit of what we can resolve numerically 
due to the extended time it takes to settle into an object 
with a well-defined maximum and minimum for the ob- 
servables. It is possible to define a 3 parameter fit to the 
data, 



(max(W)) ~ c + c ie - C2V ~\ 



while 



c + ei(l 



-c 2 Vi 



(8) 



(9) 



This parameterization can be useful in identifying the 
asymptotic value of the order parameter early in the sim- 
ulation, as we know that (0)(t — > oo) ~ Cq. It would be 



instructive, but beyond the scope of this work, to recreate 
figure 10 for a wider range of P to test if our extrapolation 
to p c indeed works as we expect it to. 

In this section we have shown examples of oscillons 
forming in various configurations from flux-tube decay 
in the 3d U(l) theory. We have also defined order pa- 
rameters to identify their presence and properties, and 
mapped out the region where we observed them numer- 
ically. Using these results, we have predicted the range 
of P where we can expect to see oscillons. For p < 0.04 
we have not found oscillons but predict their existence. 
It is a challenging numerical problem, as we have two 
very different scales in this regime (the scalar and vector 
masses) . 



V. FORMATION MECHANISMS AND FINE 
STRUCTURE 

In the previous sections we described how [7(1) os- 
cillons form directly from flux-antiflux tube interactions 
with no residual topology. On large lattices, there is a 
high probability of forming a complex flux-tube network 
which will include the needed precursor loops to form 
such oscillons. 

To better investigate the properties of these oscillons, 
even if only preliminarily, we use a simple initial condi- 
tion inspired by our 2d work on vav annihilation [19]. In 
that work, we placed a vortex and an anti-vortex nearly 
at rest and at a short distance away from each other. The 
vav pair interacts, attempts to annihilate and, depending 
on the parameters, forms a very long-lived gauged oscil- 
lon. Once this happens we put a friction wall a certain 
radius away from the oscillon to damp any outside energy 
in radiative modes that may interfere with its motion 
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FIG. 12: (Color online) 8 snapshots (t = {78, 79, 81, 83, 84, 90, 94, 106} m _1 increasing left to right) for a j3 = 0.01 oscillon. 
Plotted is a low energy electric energy- density isosurface of ~Ei ■ Ei = 0.03 (cyan) and the condensate density $t<I> = 0.5 (red). 
Each slice represents an area of 12.5 x 9mJ 2 . 




FIG. 13: (Color online) Formation dynamics of an elongated z symmetric 3d oscillon. 7i(x,y,z) — 1.6 is represented by the 
dark blue isosurface. Electric energy-density E 2 = 0.2 (cyan), and magnetic energy-density B 2 = 0.1 (green). The bottom 
configuration is llm^ 1 forward in time from the top. The plots represent ~ llm^ 1 across the longest direction. Note that 
the energy-density tube between the two oscillons in the extremities oscillates between being dominated by its electric and 
magnetic contributions. 



across the lattice. As we have seen in the previous sec- 
tions, this method allows the determination of the EFC's 
energy and structure. 

In 3d, our initial setup is similar since we can extend 
the field in the z direction without affecting the dynam- 
ics. The oscillon that forms from a flux anti-flux tube 



annihilation is just like the 2c? one except for the z sym- 
metry. To break the z symmetry we then introduce a 
spherical friction wall which effectively pinches off and 
dissipates the energy density at two polar regions away 
from the soon to be center of the 3d oscillon. Essentially, 
we are repeating the procedure of section III but imple- 
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FIG. 14: (Color online) 8 snapshots (t = {129, 143, 155, 156, 158, 164, 166, 202} mj 1 increasing left to right) for a f3 = 0.01 
oscillon. Plotted is a low energy electric energy-density isosurface of \Ei-Ei = 0.015 (cyan) and the scalar condensate amplitude 
= o.5 (red). Each slice represents an area of 25m~ 2 . The radius of the scalar condensate isosurface is approximately 
r ~ 0.75m" 1 . 



ft 




FIG. 15: (Color online) 8 snapshots (t = {133, 155, 166, 172, 178, 185, 202} mj 1 increasing left to right) for a (3 = 0.01 oscillon. 
Plotted is a low energy magnetic energy-density isosurface of ■ Bi — 0.0125 (green) and the total energy density TL = 2.5 
(blue). Each slice represents an area of 25mJ 2 . The radius of the energy density isosurface is approximately r ~ 0.75m7 . 



meriting the spherical friction wall early so as to not only 
absorb radiative modes, but to break the z symmetry 
and catalyze decay into the spherical oscillon. 

Since the stability of a time-dependent EFC is inher- 
ently linked to the exchange of energy between conjugate 
momentum (II) and potential energy (V[$]), adding fric- 
tion at the edges catalyses its decay. An example of an 
oscillonic flux tube decay catalyzed by friction is shown in 



figure 12 where a flux tube of length 12.5mJ 1 is pinched 
by a centered wall of friction of diameter d = 8m" 1 . We 
display snapshots of both the scalar-field condensate am- 
plitude (red isosurface at 4>t<fr = 0.5) and the electric 
field contributions to the energy density (light-blue iso- 
surface at \E ■ E — 0.03). Notice how the scalar field 
condensate oscillates between tube-like and bubble-like 
configurations, while the low-energy electric field alter- 
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FIG. 16: (Color online) Energy density isosurfaces for a 
/3 = 0.01 oscillon at t ~ lGOOm^ 1 after formation. Isosurfaces 
shown are at Ti = {O.Ol(blue), 0.25(magenta), 5.0(violet)}. 
The effective radius of the structure is R e ff ~ 1.15mJ , and 
the total energy in a shell of R = 4 is Z?r<4 ~ 6.9. Note 
that the majority of the energy density is nearly spherically- 
symmetric, although there is some time-dependent wobbling 
which is more pronounced at lower energy isosurfaces. 
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FIG. 17: (Color online) A few observables characterizing 
the 3d U(l) oscillon. The blue line is the maximum energy 
density of the EFC while the black line is the total energy 
within a radius of R — 4. The lower plot shows both the 

expected radius of the EFC R\ff = ^ j nff&r (s reen ) an( l 
the minimum value of $*$ within the EFC (red). Note that, 
contrary to the 2d case, the scalar field amplitude probes well 
beyond the inflection point of the double-well potential. 



nates between complex oscillatory patterns. From this 
figure, it is clear that the path to oscillon formation is 
far from trivial, even if the final configuration will dis- 
play, as we will show, near-spherical symmetry. 

It should also be noted in passing that the number of 
(nearly) independent EFC's that form is sensitive to the 
length of the initial flux tube and to where we pinch it. If, 
for example, we were to take a longer line, then, as shown 



in figure 13, nearly-spherical oscillons would form at the 
edges, connected by a flux tube between them. Perturba- 
tions during the formation process propagate down the 
structure creating an interesting pattern of oscillating 
magnetic and electric fields. Looking at the lower iso- 
surface plot in figure 13 we see that the magnetic and 
electric fields create a chain-like structure that is oscil- 
lating in both time and space. This chain then clumps 
into a certain number of oscillons per unit length. Al- 
though quite interesting, we do not pursue the study of 
these hybrid objects any further. 

Back to the oscillon structure with (3 = 0.01 of figure 
12, once the flux-antiflux line oscillon relaxes into the ap- 
proximately spherical oscillon and we have also absorbed 
much of the external radiative and thermal modes, we can 
look at the fine structure of the low energy fields to in- 
vestigate the mechanism by which energy is being slowly 
radiated at large times, as showed in figure 6. From fig- 
ure 6, we note that at large times the energy within a 
shell surrounding the oscillon can be written as 



E(t) ~ C + Cyt~ 



(10) 



where we chose a radius R — 4 around the maximum 
energy density. Extrapolating to t — *■ oo, we can see that 
the radiation will approach zero and there will be a finite 
energy left in the EFC, namely E(oo) ~ cq. Inspect- 
ing the nonspherical fine structure of the energy density, 
we find that the surviving large time structure is very 
nearly spherically symmetric and that a small fraction 
of its energy is bound up in a combination of oscillatory 
modes in the electric and magnetic fields which resemble 
transitions between excited atomic states. 

From the continuing snapshots of figure 14, we see that 
the energy density in the electric fields resembles a su- 
perposition of alternating mixed spherical harmonics. As 
time increases, the amplitude of the higher harmonics de- 
creases and the energy in the electric field becomes more 
spherically symmetric (s-orbitals dominate). If we look 
at the magnetic field contribution, as in figure 15, we ob- 
serve that the low energy component assumes a toroidal 
shape (see t = 202mJ 1 plot from figure 15). Because 
the oscillon is intrinsically time dependent, the actual 
asymptotic profile of the magnetic energy B(x) ■ B(x) is 
unclear, although we do see signs of a toroid with a con- 
tinual precession as well as an oscillation in amplitude. 
We stress, however, that these remarks are very prelimi- 
nary and that a more detailed analysis is still lacking. 

It is also important to stress that these fine structure 
harmonics only account for less than a percent of the 
actual energy of the EFC. This can be seen by considering 
that the isosurfaces in electric and magnetic field energies 
are plotted at a radius of ~ 2.5, while the effective radius 
of the structure, as calculated by the ratio 



R 



eff 



_ Jr 2 H(r)dr 
J H(r)dr 



(11) 



is much smaller (R e ff ~ 0.6 for (3 = 0.04 for example). 
In figure 16, we plot various total energy-density isosur- 



faces at late times: violet at 7i = 5; magenta at TL = 0.25; 
blue at TL — 0.01. It is clear that most of the energy is 
spherically-symmetric and concentrated in a very small 
radius, although there are departures from sphericity at 
large radii. We also present the time-dependent observ- 
ables for the 3d configuration in figure 17. This figure is 
to be compared with the 2c? version for the same parame- 
ters in figure 1 . Note how in 3d the scalar field amplitude 
probes well within the V" < part of the double- well po- 
tential. 



VI. AN INITIAL SU(2) 3d LOW p SEARCH 

In the context of non-Abelian models, the most dis- 
tinctive searches for oscillons were performed in refer- 
ences [15] and [18]. Using ansatze that approximate a 
spherically-symmetric oscillon, the authors found a sta- 
ble oscillon in the bosonic sector of the full SU(2)XU(1) 
electroweak model in 3d for the particular mass ratio 
WHiggs = 2mw- Given that this value of the Higgs mass 
is well within the reach for the upcoming LHC, the result 
has considerable relevance. Interestingly, any deviation 
from this mass ratio causes the structure to destabilize. 
We note that this mass ratio corresponds to type-2 be- 
havior in superconductor phenomenology. 

In this section, we present an initial search in the type- 
1 region where the vector particle is much more massive 
than the scalar. Of course, this is not where we expect 
to see a scalar field in the Standard Model. However, 
given that non-Abelian models are an integral part of 
any extension to the Standard Model, or of higher level 
unification, it is important to map the possible nonper- 
turbative structures that might emerge from large fluc- 
tuations about the vacuum. Only a few years back the 
Higgs itself could still have been much lighter than the 
VF-boson. Based on our U(l) work, the type-1 region is 
where we expect a coherent object to exist. Obviously, 
the non-Abelian theory is significantly different as there 
are no topological or even quasi-topological structures 
which carry as much energy as the flux tubes in the 3d 
U(l) theory. This presents a challenge to our method as 
there are no flux tubes to lock energy in the initial part 
of the simulation. 

Instead of using excessive dissipation to clean up the 
system during the initial moments of the simulation, 
we will go another route. To generate spatiotemporal 
complexity, we start the scalar fields at the unstable 
symmetric point of the potential of eq. 1, adding only 
small- amplitude thermal perturbations to excite the var- 
ious modes. As the fields evolve, they probe the un- 
stable (spinodal) portion of the potential. The ampli- 
fied instabilities - which can be viewed as isosurfaces of 
the energy density - create localized quasi-bubbles con- 
nected by quasi-strings. In Figure 18 we show a few 
observables: the red curve corresponds to the volume- 
averaged amplitude of the SU (2) Higgs field which, after 
the initial large-amplitude fluctuations, oscillates about 
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FIG. 18: (Color online) Global observables for a 51/(2) 
quench at g = 4.0 on a lattice of volume 28.8 3 mJ 3 . The 
black arrows correspond to the time-slices in figure 19. The 
system is thermalized on a Higgs potential for t = 5mJ 1 and 
then cooled with 7 = 1 [CHECK] for the same amount of 
time. After that, the system is evolved conservatively. ("I 1 ^} 
(red), (TL) (blue), and the electric and magnetic components 
of energy are plotted in cyan and green, respectively. These 
colors correspond to the isosurfaces of figure 19. 

($t<j>) r~j 0.85. The blue curve corresponds to the to- 
tal energy and the cyan and green curves at the bottom 
correspond to the "electric" and "magnetic" portions of 
the energy, respectively. It is clear that we have achieved 
good energy conservation in our simulation. The vertical 
black arrows in the figure correspond to the time snap- 
shots shown in figure 19. For clarity, the colors for the 
isosurfaces follow the conventions of figure 18. As can 
be seen in the figures, the complex network settles into 
a few large-amplitude (max [$*<!>] ~ 0) high-energy lo- 
calized objects, which then decay at about t ~ AOOmJ 1 . 
Visualizations of the simulations can be seen here . 

The fact that these structures live only to t ~ 400m -1 
says very little about their stability in near-vacuum. 
Since we have dumped so much energy into the system, 
there is a plethora of large-amplitude propagating modes 
which can and do disrupt oscillon stability. We have seen 
this in the context of U(l) models and in models with 
only real scalar fields: oscillons that have a very large 
lifetime in vacuo (t ~ lO^^mJ 1 ), if put into a high tem- 
perature background can have their lifetimes drastically 
reduced to t < 10 3 mJ 1 . This is not surprising in situ- 
ations where oscillons, which are coherent field configu- 
rations, must do work against an incoherent background 
[27]. So, it would be necessary "to clean up" the lat- 
tice which have these emergent configurations in order 
to ascertain if there is a really long lived ~ type-1 SU(2) 
oscillon. We note that in either laboratory or cosmologi- 
cal applications, symmetry-breaking is accompanied by a 
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FIG. 19: (Color online) Local observables for a SU (2) quench at g = 4.0 on a lattice of volume 28.8 3 m7 3 . Time slice isosurfaces 
corresponding to the times labeled in figure 18. Isosurfaces are twice their global averaged values, except for — 0.5. These 
particular configurations are large amplitude and decay quickly due to interactions with the propagating modes surrounding 
them. The full simulation can be viewed here . 



quench or cooling which will help oscillons - if present - 
to survive. We plan to continue this study and to investi- 
gate if there are any configurations which are long-lived 
outside of the special 2 : 1 mass ratio observed by ref. 
[15]. 



VII. SUMMARY AND OUTLOOK 

We have presented a numerical study of symmetry- 
breaking in the 3d Abelian-Higgs model. In particular, 
we searched for nontopological, time-dependent, long- 
lived configurations that may emerge dynamically as the 
system relaxes to the lower-energy asymmetric vacuum. 
Our results indicate that such oscillon-like solutions to 
the equations of motion can be easily found for a wide 
range of parameters. For the (7(1) model, the control pa- 
rameter is j3 = (m s /m v ) 2 . We found that oscillons exist 
within the type-1 regime, for (3 < 0.09. It remains to be 
seen if this range can be extended for larger (3. One short- 
coming of our approach is that it dumps too much energy 
in the initial state, as we quench from the symmetric to 
the broken symmetric vacuum. It should be possible to 



extract the approximate field behaviors that character- 
ize oscillons from our numerical solutions and use them 
as ansatze in cleaner searches as, say, was done in refs. 
[14] or [15]. This way, we could probe into the type-2 
regime in search of oscillons. Our preliminary study of 
their rich resonances and of possible oscillons in type-1 
5/7(2) models suggest that there is much more to be ex- 
plored. Although analytical results are very challenging, 
it is possible that some may be achievable extending the 
procedure of ref. [24] to models with gauge fields. 

The fact that oscillon-like configurations emerge spon- 
taneously during symmetry breaking should not be over- 
looked. There is very poor understanding of the dynam- 
ical aspects of symmetry breaking in gauge theories as 
they involve complex nonlinear and possibly nonequilib- 
rium physics. The emergence of oscillon-like EFCs in the 
Abelian-Higgs model indicates that the thermodynamics 
of these systems is far from trivial. In cosmology or in 
colliders, the presence of these configurations will delay 
equipartition and thus the final approach to equilibrium. 
They may, for example, affect the calculation of the re- 
heating temperature in inflationary models [10]. They 
may also affect the structure of the vacuum and thus 
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the computation of transition amplitudes between vac- 
uum states [1]. We are now entering an era where even 
desktop computers have enough power to perform such 
simulations. It is clear that much new physics remains 
to be explored. 
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APPENDIX A: LATTICE IMPLEMENTATION 

In order to construct gauge invariant objects in a 
Hamiltonian lattice formulation we will be using the Wil- 
son loop formalism [28] . For the lattice Laplacian, we use 
the second-order version of the covariant derivative [29] 



r/tM <f>- „ 



> x }+(D(6x 2 ), 
(Al) 

where the hatted variables x denote spacctimc lattice po- 
sitions and p denote a lattice displacement by one lattice 
spacing 5x. The gauge field strengths are constructed 
from the unitary link variables U(x, p) satisfying the lat- 
tice gauge transformations 



U{x,n) = Ug 



A-U^A^ 



(A2) 



We can use the link variables to construct the unitary 
tensor 



U l Mx) = IPuLuiJu^, 



(A3) 



which is the fundamental lattice gauge invariant object 
from which to derive the field strengths. Here, i,j repre- 
sent spatial indices. The magnetic contributions to the 
Hamiltonian can be found from the combination 



1 



F 13 ■ Fa, 



1 



1 - Re 



Tr(C/g") 



, (A4) 



The same result can be obtained from constructing the 
covariant derivative of the field strength in the adjoint 
representation by taking a gauge invariant backward 
derivative on the field strength tensor, 



g5x 3 



Im 



{u^(x)-U^U^(x-0)U^}, 

(A8) 

where for an Abelian theory the extra U~_^s drop out 
to unity. 

The momentum components are just simple complex 
matrices at each point in space for the scalar, and a vec- 
tor of complex nonunitary matrices for the electric fields. 
Since we are in a gauge in which the temporal vector 
potential is set to zero, derivatives with respect to time 
are simply, Ei — d t Ai. We can then write the gauge 
equations of motion on the lattice as, 



d t El 
dtAi 



3 _ 



EL 



Ji 



where the current is 



2£ 

8x 



Im 



( < i , :.'.:>-, 



(A9) 
(A10) 



(All) 



This equation is now in a perfect form for our higher- 
order symplectic integration scheme, and works also with 
the link formalism, provided that we change equation 
A10 to the form, 



W(x,t+) = Wix^e- 195 * 5 ^ 3 ^'^. 



(A12) 



That is, provided our numerical exponential is accurate 
(and the electric field stays within its group during simu- 
lations with stochastic forcing, see below) then the links 
stay unitary at all times. The gauge condition A n = 
giving the time component of equation 3 becomes a non- 
dynamical constraint equation, the Gaussian constraint, 



^(4-4^=2^(4$). 



5x 



(A13) 



where p = 2 for SU(2) and p 
strength is 



1 for U(l). The field 



p l 3 



1 



T m TJ % 3 



(A5) 



The scalar field will be a nonunitary complex matrix 
which transforms as 



$(Z) EE -> A £ <f>£ 



(A6) 



The covariant derivative of the force tensor needed for 
the equations of motion can be found from the variation, 



SAJy) 



(A7) 



Provided that this is satisfied initially, then it is main- 
tained to numerical accuracy at all times for the U{V) 
theory. There are some issues in exact convergence for 
the non-Abelian theories. 

The scalar equations are also split into conjugate mo- 
mentum fields to give 

d t n x = (Ai4) 

d t f> x - Ha (A15) 

where T)f att $ is defined in eq. Al. The Hamiltonian 
density Tl x that corresponds to these equations is then 
given by, 



19 



u, = nj, • + vM ■ vvi>., + ' /•;, • e* 



q 2 Sx 4 ' 



) + 



i (V. • 



(A16) 



The Langevin implementation of these equations is ex- 
plained in the next two appendices. 



APPENDIX B: GAUGE 
LANGEVIN IMPLEMENTATION 1: FORCING 
MOMENTA 

In this appendix we summarize the two methods of 
implementing stochastic noise and friction terms in sys- 
tems with local gauge symmetries and thus Gaussian con- 
straints. The methods differ in that whereas one couples 
noise to the fields, the other couples noise to the conju- 
gate field momenta. We also explain how the implemen- 
tation which forces the fields can be made to satisfy the 
Gaussian constraint equation at all times. Conversely, 
the method of forcing the momenta destroys this con- 
straint and thus the system must be "stopped," (that is, 
one must set all momenta set to zero) before the true 
simulation (symmetry breaking) begins. 

It is useful at times to run simulations with only the 
dissipative terms present and no thermal noise in order to 
relax to a ground state or specific configuration. We also 
present a proof which shows that spatially homogeneous 
friction does not affect the Gaussian constraint and is 
safe to use in a simulation. 

Consider a U(l) gauge field coupled to a complex scalar 
field with some Hamiltonian H = Ti. which is a func- 
tion of the conjugate momentum fields for the scalar 

and the vector, 7r Q = |n, e\, and the respective fields 

(j) a = |<i>, A\. The Langevin dynamics for this system 
is, 

dtK + iK = -a*.l>«,0S) + a- (bi) 

X 



where £ = y ^^ rndcO, where rndc() is a Gaussian 
random number of zero mean and unit variance. Since 
7T° = II is a complex number, we must add a random kick 
for each degree of freedom and thus £ v o = £p^ e + iClm i 
where both the real and imaginary components of are 
identical except for their random numbers. Numerically, 
they are just sequential calls to the random number gen- 
erator. If we had written the set of equations as a scalar 
field of 0(n) symmetry then we would just call the ran- 
dom number £ n times and couple that force linearly to 
each component of the 0(n) field. This is also the same 



for the vector field as each ir 1 = E l component gets its 
own random kick. 

This Langevin evolution violates locally the Gauss con- 
straint, 

i J2(K-i - K) = i(*°M? - ("SM)- (B2) 

i 

Since this equation is preserved by the equations of mo- 
tion without the forcing noise, to satisfy it we can just 
set all the momenta to zero at the beginning of a simula- 
tion. Although this procedure also destroys the thermal 
state, it will approximately be restored after the energy 
is quickly transfered from the spatial derivatives to the 
temporal ones as the system tries to reach equipartition. 
We stress that we are not interested in enforcing a strict 
thermal state anyway, only in using the Langevin ap- 
proach to excite a large spectrum of field modes. 

Note that the 7 frictional coefficient is not a problem 
as it scales smoothly in &-space as e -7 ' and so does not 
change any spatially local term. We can show that this 
has no effect by considering the equations 

<9 f vr + 77r = V 2 <t)-m 2 eff (j) (B3) 
d t E l +~/E l = diF^ + f (B4) 
d t E l = J°. (B5) 

Taking the divergence of the second equation and the 
time derivative on the last and then subtracting the two 
gives: 

+ jd t E l = O + diX-dtJ . (B6) 

In order to get the appropriate correction that should 
cancel the 7 term we expand 

d t J° -» 2glmd t ^4> (B7) 
-> 2 fl Im((-77r t +V 2 ^ - m 2 eff tf)(f>) (B8) 

-► 0^-25710171-^. (B9) 

So equation B6 can be expressed as 

idiE? = 7 J°, (BIO) 

which means that the continuity equation is of the same 
form as the initial 7 = constraint equation B5. 

APPENDIX C: GAUGE 
LANGEVIN IMPLEMENTATION 2: FORCING 
FIELDS 

Recently we became aware of the work by Krasnitz 
on how to maintain constraints on gauge fields during 
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Langevin evolution [30]. After some experimentation 
with the algorithms in reference [30] , we were able to find 
a simpler and less computationally demanding method of 
maintaining the Gaussian constraint on the ?7(1) gauged 
scalar system. While not as general as the methods of 
[30], for at least a subset of the Abelian systems our 
method is simpler and numerically faster Langevin-type 
thermalization method. In Krasnitz's method, the force 
function must be called twice per time step and the sim- 
ulation time almost doubles. 

For clarity, we repeat the basic equations for an 



Abelian system with friction, 

d t A t = Ei (CI) 

d t $ = n (C2) 

d t Ei = djF ji --fE i + 2gIm<S>V i <S>* (C3) 

d t Il = X> 2 $-7lI- V'[$] (C4) 

and the Gaussian constrain equation 

d i E i = 2glmn*$., (C5) 



where it is implied that we use the lattice version of the 
operators. 

The method works as follows: instead of adding a 
stochastic force which satisfies a fluctuation-dissipation 
relation to the momentum evolution equations, we add 
a random perturbation to the field equations. Thus, the 
amplitude perturbations are written as 

8 t Ai = E. + Ce (C6) 

d t $ = n + Cn. (C7) 

It is obvious that a change in amplitude of A{ cannot 
affect the electric field and so the divergence in equation 
C5 is not modified. Let us look at the charge density 
at two separate times to see what conditions on £n will 
maintain the right side of the constraint equation C5. 
Writing $(< + St) = $+ = $ + USt + ( n St, gives the 
charge density at half a leapfrog step, 

J°($+) = J°($ + ILSt + tnSt) (C8) 
2 ff Imir$ + = 2gImIL*($ + n5t + (nSt) (C9) 
= J°+2gMmn*Cn, (CIO) 

where we only need to advance <& — > because of the 
leapfrog scheme of integration. If J is locally conserved 
at this half step then although it is modified in the next 
half step, it is not changed in any way which violates 
equation C5. Withis we can get the necessary form of (ji 
which makes ImlFCn = and thus maintains the local 
Gauss constraint. If we choose a real stochastic variable 
C and define 

c " s V5FH C > (C11) 



we can add random fluctuations at any time during the 
simulation without having to set all momenta to zero to 
regain the Gaussian constraint as we must with the first 
type (momentum forcing) of implementation. 

It should be noted that this particular implementation 
will not randomize the phase of an ungauged complex 
scalar, since it ends up being only a density perturbation 
which relics on the gauge field to randomize the phase. So 
to thermalize a complex scalar, use a complex £ = 0?+*Cf 
instead of Cn ■ 



Criticisms and Challenges 

There is one obvious criticism to the thermalization 
technique presented in this section. If we look at the 
equations, we see that we are in effect only modifying 
the spatially-dependent degrees of freedom by making 
fluctuating charge and magnetic densities but not do- 
ing anything that would make the charge-anticharge and 
electric distributions come to equilibrium. 

For example, in reference [30] the forcing does change 
both sides of equation B2 but in a way that is balanced. 
This seems a better way of doing things at the cost of 
a few more computational steps, as it does not change 
the equilibrium configuration that the forcing directs the 
system towards. Just because we are not modifying ei- 
ther side of equation B2 does not mean that the charge 
densities remain trivial. Energy flows very quickly and ef- 
fectively from the current to charge-anticharge and from 
magnetic to electric degrees of freedom inherently be- 
cause of the equations of motion. So, a thermal equi- 
librium distribution of diEi and spatially varying charge 
densities does in fact form. 

Another issue to be careful of is that when thermal- 
izing from the amplitude perturbations, the largest one 
can make the electric contributions to the temperature 
is (EiE 1 ) = (BiB 1 ), which comes from minimization of 
the Lagrangian instead of distributing % per degree of 
freedom in the Hamiltonian. 

This happens anyway if we maintain the Gauss con- 
straint by setting all the momenta to zero. None of these 
techniques give a perfect thermal distribution. 

With both techniques we see that generally the mo- 
menta thermalize well (11*11) oc {EiE 1 ) oc T as do the 
surface kinetic terms (di$*d l &) . But the energy mixture 
can have issues between the magnetic and the electric 
terms. These issues are present in every form of real 
time classical thermalization that we have tried. 
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